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EVALUATION 


The  present  contract  provides  basic  theoretical  and  experimental  in¬ 
formation  on  some  of  the  fundamental  material  properties  of  CCD  de¬ 
vices  such  as  noise  sources  due  to  surface  states,  and  high  field  and 
interband  ionization  effects  which  limit  the  performance  of  such  de¬ 
vices  in  implementing  signal  processing  functions.  The  results  ob¬ 
tained  will  be  useful  in  enhancing  the  performance  of  and  defining  the 
limitations  of  such  devices  as  to  speed  of  operation,  dynamic  range, 
and  functional  density. 

Related  TPO  //  and  Title:  A  -  Technology/D  -  Solid  State  Devices 
Sub-Thrust  //  and  Title:  2  -  Signal  Processing  Devices _ 
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I .  INTRODUCTION  1 

A  number  of  fundamental  material,  properties  which  limit  the 
performance  of  silicon  charge  coupled  ievices  for  signal  processing 
are  investigated  in  this  program.  To  get  a  wide  dynamic  range,  the 
noise  sources  must  be  identified  and  reduced.  To  increase  the  frequency 
of  operation,  the  basic  limitations  on  the  charge  transfer  speed  and  charge 
transfer  efficiency  must  be  assessed.  This  report  gives  a  summary  of 
the  studies  which  aimed  at  understanding  some  of  these  basic  limitations. 

To  facility  the  study  of  noise,  especially  at  low  frequencies  in 
MOS  transistor  and  CCD  structures,  a  hardware  computer-controlled  noise 
correlation  processor  was  built.  Its  design,  error  analysis  and  operation 
were  described  in  a  published  RADC  report  [1]  and  a  summary  is  give  in 
section  II  in  this  final  report. 

The  1/f  low-frequency  noise  is  generated  mainly  in  the  output  MOS 
transistor  stage  of  a  CCD.  For  the  CCD  transfer  stages,  the  dominant 
noise  source  is  the  random  thermal  generation  and  capture  or  trapping  of 
electrons  (in  n-channel  CCD)  at  the  SiO^-Si  interface  states.  A  detailed 
review  of  the  CCD  noise  literature  is  given  in  a  second  RADC  report  [2].  In 
section  IV,  a  summary  of  the  interface  trapping  noise  is  given  which  extends 
previous  work  on  thermal  emission  noise  to  random  capture  noise  as 
well  as  noise  due  to  incomplete  filling  of  the  interface  states  for 
low  and  moderate  capture  rates. 

The  ultimate  transfer  speed  in  CCD  is  limited  by  two  dominant  high 
electric  field  effects  which  occur  when  the  gap  and  gate  lengths  are 
reduced  to  shorten  the  charge  transfer  time.  These  are  the  mobility 
reduction  due  to  the  increasing  electron  energy  loss  to  phonons  in  high 
electric  fields  and  the  interband  impact  ionization  or  electron-hole  pair 
production  in  high  electric  fields  which  obscures  the  zeros  and  which 
also  contributes  noise. 


2 


In  section  III,  the  effect  of  high  field  mo  b  i '  It .  reduction  on 
the  charge  transfer  speed  ’  s  ir.^I/zed  and  the  use  of  appropriate  hot 
electron  ^ltfusivity  is  discussed.  The  ultimate  transfer  speed 
limitation  by  high  electric  fields  is  discussed  in  section  VII  which 
shows  its  dependence  on  gate  length,  mobility  reduction,  drift  velocity 
saturation  and  impact  ionization  effects. 

A  third  high  field  effect  is  actually  beneficial.  This  is  the 
reduction  of  the  thermal  capture  rate  of  electrons  at  the  interface 
states  in  high  fields.  This  would  reduce  the  interface  trapping  noise 
as  well  as  transfer  efficiency.  It  is  discussed  in  section  IV. 

Fundamental  measurement  methods  are  described  in  section  V  and  VI 
to  characterize  the  interface  states  in  surface  channel  CCD's  and  bulk 
deep  levels  in  buried  channel  CCD's.  The  small-signal  equivalent  circuit 
model  is  described  in  section  V  and  is  used  to  obtain  the  density  of 
states  and  the  thermal  emission  and  capture  rates  of  electrons  and  holes 
at  the  interface  states  from  MOS  admittance  measurements.  In  section  VI, 
the  voltage  stimulated  capacitance  transient  spectroscopy  (VSCTS),  an 
improved  DLTS  (Deep  Level  Transient  Spectroscopy),  is  described  as  a 
method  to  determine  the  recombination  parameters  of  electrons  and  holes 
at  impurity  levels  and  ion-implantation-damage  levels  in  the  bulk  of 


silicon. 


II.  NOISE  CORRELATION  PROCESSOR 

2. 1  Introduction 

A  hardwired  noise  correlation  processor  has  been  designed,  built 
and  tested.  The  correlator  consists  of  a  high  speed  noise  data 
sampler  and  multiplier.  Noise  data  can  be  sampled  at  selected  sample 
rates  (up  to  200  KHz)  and  selected  number  of  samples  (up  to  16384 
samples).  These  data  are  then  processed  by  the  high  speed  multiplier 
and  stored  in  the  memory  of  the  processor.  They  are  then  clocked  into 
the  memory  of  the  Wang  720C  programmable  desk  calculator  for  further 
analysis  and  plotting.  The  details  of  design,  error  analysis  and 
operation  of  the  noise  correlator  was  reported  in  a  previous 
report  [1].  The  results  are  summarized  in  the  next  sub-section. 

2.2  Experiment  Results 

A  detailed  block  diagram  of  the  correlator  is  shown  in  Figure  2.1. 
The  instrument  has  several  convenient  features  such  an  input  overload 
and  output  ready  indicators,  adjustable  clock,  adjustable  output  gain, 
adjustable  number  of  samples,  and  an  output  recirculate  mode  for 
oscilloscope  display.  Oscilloscope  photographs  of  actual  correlator 
calculated  autocorrelation  functions  ure  shown  in  Figure  2.2  for  sinu- 
sodal,  square  wave,  white  noise  and  Lorentzian  noise  inputs.  The  upper 
trace  in  each  case  is  the  autocorrelation  function  and  the  lower  trace 
is  the  input  signal. 

Figures  2-3  and  2-4  show  plots  of  the  correlator-calculated  auto¬ 
correlation  function  and  the  subsequent  calculator-calculated  power 

2 

spectrum.  The  input  signal  here  is  Lorentzian  noise  with  a  l/[l+(f/f  )  ] 


spectrum.  There  is  some  scatter,  but  a  three  point  average  shown  in 
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Figure  2.5  gives  a  smooch  curve  with  liccle  loss  in  frequency  resolution. 
In  fact,  Che  small  filcer  hump  ac  about  600Hz  is  still  retained. 


2 . 3  Summary 

Various  sources  of  error  were  considered  in  designing  the 
correlator  in  order  to  get  the  most  efficient  system  in  terms  of  speed, 
accuracy,  size,  cost  and  convenience.  These  included  high  frequency 
aliasing,  low  frequency  resolution,  A/D  conversion,  finite  time, 
finite  sampling  and  quantization.  The  final  correlator  uses  the  incomplete 
correlation  method  and  four  bit  input  resolution.  This  system  combines 
the  high-speed,  dedicated  front-end  processor  (the  correlator)  with  a 
flexible  general  purpose  computer  which  can  be  a  calculator,  minicomputer 


or  microcomputer. 


Out 


(r)  (J) 


Figure  2.2  Oscilloscope  photographs  of  the  correlator  computed 
autocorrelation  functions  for  (a)  sinusoidal,  (b) 
square  wave,  (<•)  white  noise  and  (d)  l.orentzian 
noise  input  signals.  The  input  signals  are  shown 
in  the  lower  part  of  each  figure  and  the  autocorrelation 
functions  are  shown  in  the  upper  part  of  each  figure. 
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Figure  2.3  The  autocorrelation  function  of  a  Lorentzian  noise 
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2.4  The  power  spectrum  of  the  Lorentzian  noise  shown  in 
Figure  2.3. 


III.  IOT  ELECTRON  EFFECTS  IN  CHARGE  TRANSFER 

3.1  3robIem  Identification 

The  charge  transfer  process  in  short-gate  surface-channel 

charge- coupled  devices  is  mainly  governed  by  the  fringing  electric 

4 

field  wnich  can  be  as  high  as  5  x  10  v/cai  [3,4].  Electric  fields  of 
this  magnitude  heat  up  the  charge  carriers  to  kinetic  energies  of 
several  thousand  degrees,  far  above  the  lattice  temperature.  This  is 
known  as  the  hot  electron  effect  (5].  One  of  the  consequences  of  the  hot 
electron  or  high  field  effects  is  the  reduction  of  mobility  and 
diffusivity  due  to  the  increased  scattering  rate.  This  could  cause  drift 
velocity  saturation  and  limit  the  charge  transfer  rate.  Trapping  noise 
due  to  the  interface  states  in  CCD  devices  is  also  affected  by  the  presence 
of  the  high  electric  field.  The  high  electron  temperature  or  the  shifted 
electron  distribution  in  energy  in  the  conduction  band  will  drastically 
affect  and  usually  reduce  the  thermal  capture  rate  of  electrons  at  the 
interface  states.  The  high  electric  field  will  also  increase  the  thermal 
emission  rate  of  trapped  electrons  at  the  interface  states  due  to 
lowering  of  the  trap  potential  barrier  by  the  electric  field. 

In  addition,  due  to  the  high  carrier  temperature  and  electric  field, 
it  is  possible  to  cause  interband  impact  generation  of  electron-hole 
pairs.  This  can  impose  limits  on  the  gate  and  gap  lengths  and  gate 
voltages  of  CCD  structures.  Trapping  noise  is  considered  in 
section  IV.  A  computer  model  was  developed  to  provide  a  quantitative 
evaluation  of  the  hot  electron  effect.  This  is  described  next. 
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Figure  3.1  The  CCD  structure  used  in  the  computer  model. 


3.2  Theoretical  Model 

The  structure  of  the  transfer  gates  used  in  the  model  is  shown 
in  Figure  3.1.  We  are  only  concerned  with  charge  transfer  between 
transfer  gates.  The  input  and  output  stages  of  the  CCD  are  not 
considered  here.  The  relevant  equations  are: 

(i)  The  continuity  equation 


3n  3J 

q  IT  =  TT  • 


(3.1) 


(ii)  The  current  equation 


J 


+  qpnE. 


(3.2) 


(iii)  The  Poisson  equation 
2 

H  »  0,  for  0<  z<  z  (3.3a) 

—  -  ox 

2 

£g7  i>  =  q  (n  -  Nj.),  for  z  <  z  <  »  (3.3b) 

Equation  (3.13  through  (3.3)  are  solved  simultaneously  subject 
to  the  boundary  conditions  given  in  reference  [3].  Equation  (3.2)  is 
uporcpriate  for  describing  hot  electron  phenomena  and  was  first 


derived  by  Stratton  [6],  who  assumed  that  the  scattering  rate  does  not 

depend  on  the  position  x.  In  the  case  of  a  position  dependent  scattering 

rate,  the  approach  of  Sah  and  Lindholm  has  to  be  used  [7]. 

The  parameters  in  equations  (3.1)  through  (3.3)  are  normalized 

according  to  Table  3-1.  This  is  equivalent  to  setting  e  =  1  and  q  =  1 

s 

in  Equation  (3.1)  to  (3.3).  In  Table  3-1,  Uq  and  are  the  zero  field 
mobility  and  diffusivity.  The  numerical  values  of  the  parameter  used 
in  the  calculations  are  also  given  in  Table  3.1. 


Table  3.1  Definition  of  Normalized  Symbols 


Parameter 

Symbols 

Normalized  to 

Value 

Length 

X,  z 

>  k„T, /e^N 

D  s  B  L  A 

1.30E+3  A 

Time 

t 

XDe/,;jokDTL 

9.34E-12  s 

Potential 

4 

kBTL//e 

2.58E-2  V 

Concentrations 

n.p.Nj, 

Vd 

1.30E+10  cm-3 

Electron  Temp. 

T 

tl 

3.00E+2  K 

Mobility 

u 

u 

0 

7.00E+2  cm2/V-s 

Diffusivity 

D 

D 

o 

1.81E+1  cm2/s 

To  further 

simplify  the  calculation,  we 

approximate  the 

inversion  layer  by 

a  sheet  of 

electrons.  Thus, 

Equations  (3.3a) 

and  (3.3b)  become 

o 

II 

O 

CM 

t> 

0  <  z  <  z 

ox 

(3.4a) 

7^-$  =  qn(x)5(z) 

z  <  z  <  b 

ox 

(3.4b) 

*  -  0 

b  <  z  <  00 

(3.4c) 

Here,  b  =  z.  +  z 

d  ox 

and  z,(t) 
a 

is  the  thickness 

of  the  silicon  surface 

space  charge  layer 

under  the 

gate  where  the  charges  are  transferred. 

The  advantage  of  Equation  (3.4)  over  (3.3)  is  that  (3.4)  can  be 
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solved  using  Fourier  expansion  of  $  instead  of  the  finite  difference 

method  so  that  substantial  saving  in  computing  time  is  achieved.  Our 

results  for  the  fringe  field  have  been  compared  with  the  exact  results 

of  Carnes,  Kosonocky  and  Ramberg  [8].  It  is  found  that  the  minimum 

fringe  field  as  obtained  from  Equation  (3.4)  agrees  within  20  percent 

with  the  exact  results  from  using  equation  (3.3).  Elsewhere,  the 

fringe  field  values  are  closer  to  the  exact  solution  than  the  minimum 

value  case.  Two  additional  simplifications  are  used.  They  are  described 

as  follows.  The  potential  between  the  gates  at  z  =  0  was  assumed  to 

vary  linearly.  Therefore,  potential  barriers  do  not  exist  between  the 

gates.  This  approximation  is  good  for  gap  length  <lum  [9]  as  considered 

here.  For  larger  d,  the  approximation  generally  breaks  down.  The 

second  simplification  is  that  z^ft)  is  taken  to  be  a  constant  parameter. 

The  value  for  z^ft)  does  not  affect  the  results  significantly  for 

short-gate  CCD's.  The  results  presented  change  less  than  2  percent  if 

z  is  varied  by  a  factor  of  two.  The  actual  values  of  z,  is  determined 
d  a 

by  and  n(x, t) . 

In  order  to  obtain  the  numerical  solutions  of  the  system  of 

equations  given  by  (3.1),  (3.2)  and  (3.4),  the  electric  field 

dependencies  of  the  mobility  and  diffusivity  must  be  known.  The 

diffusivity  is  given  by  D  =  (kgTe/q)u,  where  Tg  =  Tg(E)  is  the  electron 

temperature  of  the  Maxwellian  distribution  and  u  =  u(E).  The  diffusivity 

does  not  depend  strongly  on  electric  field  since  a  decrease  of  mobility 

with  electric  field  is  compensated  for  by  an  increase  of  electron 

temperature.  Thus  in  our  calculations,  we  take  D(E)  =  D.,,  t.he  zero 

u 

electric  field  diffusivity.  This  is  a  good  approximation  for  the 
diffusivity  value  measured  in  the  bulk  [6,  10,  11].  Elaid  et  al  [4] 
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used  a  diffusivicy  given  by  D  =  (k  T  /q)u  where  T  is  the  lattice 

B  L  L 

temperature.  This  clearly  underestimated  the  diffusivity  when 

T  >>  T  . 
e  I, 

The  electric  field  dependence  of  mobility  can  be  calculated  from 
che  following  power  balance  equation  where  e  equals  electron  kinetic  ener 

<Je>  »  <3e/3t>n  (3.5) 

The  rate  of  energy  loss  to  the  lattice  <3c/3t>  for  low  electric 
field  Tg  =  T^  and  for  T^  =  300  K  may  be  approximated  by  the  two- 
dimensional  formula,  since  we  approximated  the  electron  distribution 
in  the  silicon  surface  space  charge  layer  by  a  two-dimensional  sheet 
as  in  Equation  (3.4b).  This  equation  was  given  by  [12] 

<3e/3t>  =  2 . 8xl0~9  [(Tg/TL  -  1]  Watt  (3.6) 

where  the  bulk  electron- intervalley  phonon  coupling  constant  is  used. 

In  the  region  of  high  electric  field  when  T  >  3  T  ,  an  alternate 

£  -L 

formula  is  used.  The  high  field  formula  is 

<3e/3t>  *  9 . 64xl0~9  /Tr  Watt  (3.7) 

e  L 

This  is  a  three-dimension  formula  for  energy  loss  to  intervallev 
phonons.  The  electron- intervalley  phonon  coupling  constant  from 
Jacoboni  et  al  [11]  was  used.  The  spatial  dependencies  of  electron 
temperature,  T  ,  are  computed  from  Equation  (3.5),  (3.6),  or  (3.7). 

In  addition,  we  need  the  dependence  of  mobility  on  electron 
temperature  in  order  to  obtain  its  electric  field  dependence.  The 
explicit  formula  is  a  complicated  integral  over  the  Maxwellian 

distribution,  exp(-E/k  T  ),  where  E  is  the  kinetic  energv  of  t.e 

Be 


electron,  since  both  the  acoustic  and  intervallv  phonon  scattering 
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are  important  in  determining  the  momentum  or  mobility  relaxation  time. 
This  is  further  complicated  by  scattering  from  impurities  and  oxide- 
silicon  interface  imperfections.  To  use  the  explicit  formula  wouli 
increase  the  computation  time  by  one  order  of  magnitude.  Instead  .n.<- 
following  approximation  is  used  for  u(E) 

u  -  uQ  [1  +  (E/Ec)V 1/2  (3'8) 


Here  E  is  the  critical  field.  This  is  a  good  approximation  at  low 
c 

field  when  E  <  E  and  the  electric  field  effects  on  charge  transfer 
c 

and  noise  are  not  important.  It  also  predicts  the  observed  velocity 
saturation  at  high  electric  field,  wl  en  E>>E^_. 

The  critical  field  was  determined  by  numerical  estimates  which 
showed  that  ^  at  E  =  E£  is  about  2T^.  Thus  taking  D  =  DQ,  Tg  =  2Tl> 
and  using  Equations  (3.2),  (3.5),  (3.6)  and  (3.8),  we  find 
E  =  -  (sign  E) • (k„T  / /2qn) * (  in/ 9x) 

+  k(kgT3n/9x)  2  /2qn  +  2  .  HxlO-9  /  (qMg/*^)'  (3.9) 

For  very  high  surface  potential  values,  E^  also  depends  on  the  surface 
potential  [14].  This  dependence  is  not  considered  here.  Energy- 
transport  in  the  x-direction  is  neglected.  This  is  important  oni  ■ 
for  extremely  large  values  of  3T(a/3x. 

The  system  of  equations  was  solved  by  the  following  iteration 
procedure.  The  Laplace  solution  of  the  electrostatic  problem  is 
first  obtained  by  setting  n(x)  =  0  in  Equation  (3.4),  and  using  the 
Fourier  series  expansion  of  the  potential  given  by 


:: (x.z.t) 


■(t)*exp(tk  c)*sin(k  x) 


(3.10a) 
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This  is  similar  Co  Che  mechod  of  reference  [8]  buC  ic  differs  in  chac 

they  considered  only  che  case  of  z  —  ».  The  solucion  of  Che  Poisson 
9  1 

equacion,  n  ^  =  (q/e  ) n (x , c )5 (z )  is  added,  which  is  approximaced  by 

[3,  14] 


5l(x,Z,t)  -  -qn(x,c) 


(z  U  )  (zA  ,) 

ox  ox  d  si 


(z  /f 
ox 


(z./z.) 
a  si 


(3.10b; 


=  -qn(x,c)*(z  /,;  )  (3.10c) 

ox  ox 

The  approximaCion  given  by  Equacion  (3.10c)  holds  if  (z  It  )<<(z,/c  . 

°  ox  ox  d  si 

This  3pproximacion  reduces  che  Cransfer  time  by  5  percenC  ac  mosc  buc 
ic  considerably  reduces  che  conpucacion  cime  since  Che  boundary  z  =  h 
used  in  che  charge-free  solucion,  sJ(x,z,c),  is  then  independent 
of  ^(x.z.c).  The  pocencial  is  given  by 

+>(x,z,t)  =  ?5(x,z,c)  +  (p'(x,z,c)  (3.11) 

carrying  ouc  Che  calculations,  it  was  necessary  to  take  more 
than  400  terms  in  Che  Fourier  series  expansion,  Equation  (3.10a),  to  obtain 
0.1  percent  accuracy  in  the  potential,  $ .  The  electric  field  is  then 
calculated  from  E  =  -yj/qx  using  the  solution  given  by  Equation  (3.11). 

This  is  used  in  Equations  (3.1)  and  (3.2)  which  are  then  solved  by  the 

finite  difference  mechod  as  described  in  [4j.  The  electron  density 
at  the  next  time  step,  n(x,t  +  £t)  is  then  obtained.  It  is  used  in 
Equation  (3.10c)  to  give  a  new  (x,z, t+^C) ,  and  a  new  total  potential, 

:  i  ,z,C-  t).  The  new  potential  is  then  used  to  calculate  the  electric 
field  and  electron  density  for  the  next  time  step.  This  procedure  is 
repeated  until  all  the  charges  in  the  initial  distribution  n(x,t=0)  are 
transferred  from  one  gate  to  the  adjacent  gate.  This  requires  about 
3  x  IQ4  normalized  time  steps. 
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In  several  calculations ,  the  final  steady-state  carrier 
distribution  or  potential,  > ^ (x, z , t=co) ,  was  used  as  the  initial 
carrier  distribution.  However,  an  approximate  initial  carrier 
distribution  was  found  to  give  better  that  5  percent  accuracy  in 
the  computed  transfer  time  for  the  smallest  gap  and  gate  used  which 
are  0.5  and  1  cm  respectively.  The  accuracy  is  better  for  the  other 
parameters  and  or  the  transfer  times  of  larger  gap  and  gate  lengths. 
This  approximate  intiai  charge  distribution  has  an  edge  distribution 
of  n(x)  =  2n(0)/[l  +  exp(x/A)]  where  A  =  (L+d)/30  and  x  =  0  is  the 
edge  of  the  gate.  Most  of  the  numberical  computations  made  use  of 
this  initial  charge  distribution  in  order  to  reduce  the  computation 
time  by  a  factor  of  two.  Finally,  it  should  be  pointed  out  that  it 
is  incorrect  to  use  y  3  u(E)  for  all  times  as  was  done  in  [4].  As  f*  ”  , 
Equation  (3.9)  gives  a  very  large  E^  so  that  Equation  (3.8)  gives 
y  =  u~  at  the  end  of  the  transfer  when  J  =  0  and  T  1  T  . 

U  6  Li 

3-3  Numerical  Results  and  Discussion 

The  following  values  for  the  various  device  parameters  were 
used  in  the  numerical  calculations:  zqx  =  10-5  cm,  d  =  0.5xl0~4  cm. 

TL  =  300  K’  ESi  =  U'7e0  and  £ox  “  3 ’ 9e0’  Zd  =  b  ~  Zox  ”  10  4  Cm> 

2  o 

Uq  =  700  cm  /v-s,  Dq  =  18.1  cm’/s.  The  saturation  velocity  was 
b 

■*  x  10  cm/s  and  the  difference  in  gate  voltage  between  adjacent 
electrodes  during  transfer  is  10  V. 

Figure  3-2  shows  the  transferred  charge  densitv  as  a  function 
of  position  along  the  interface  at  four  times.  Figure  3.3  shows  the 
charge  density  as  a  function  of  time  at  three  positions.  Figure  3.^* 
shows  the  coefficient  n  versus  transfer  time.  h  is  defined  as  the  ratio 
of  charge  left  behind  under  the  storage  gate  to  the  initial  charge 


Fig.  3.2  Transferred  charge  density 
as  a  function  of  position 
along  the  interface  x. 
t^=0. 0025ns,  t2=0. 0125ns, 

t  =0. 0275ns,  t  =0. 0775ns. 

3  4 


Fig.  3.3  Transferred  charge  density 
as  a  function  of  time. 
X2*-0.75ym,  X2=0.00pm 

x3=+0.75pm 


Fig.  3.4  Fractional  residual  charge 
under  the  gate  for  n=2.5xlO 
cm-2  . 

Curve  1  y  =  y„ 

D  =  D0 

Curve  2  y  =  y(E) 

D  =  (kBTL/q)y(E) 
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under  Che  same  gate  when  no  Interface  states  are  present.  Curve  1 

neglects  the  heating  of  the  charge  carriers  completely.  Curve  1  is 

calculated  after  using  the  formulation  of  reference  [4]  where  hot 

carrier  diffusion  is  not  included  properly,  and  Curve  3  is  calculated 

as  described  in  the  previous  section  using  a  constant  diffusivity 

D  =  D(E=0)  =  Dq.  The  largest  effects  are  caused  by  the  mobility 

reduction  as  already  noted  in  reference  [4].  Figure  3.5  shows  the 

-4 

electron  temperature  under  two  transfer  gates  of  10  cm. 

I  GATE _  ■  lO  ■J  _ GATE _ j 


Figure  3.5  Electron  temperature  below  the  transfer  gate  for 
a  gate  length  of  1  urn  at  two  transfer  times. 
Carrier  concentration  is  n»2 . 5x10 1 1  cm-2 . 

At  thermal  equilibrium,  the  mass  action  law  requires  that 

=  c^  N  exp(-  /k_T) .  3ut  e°  =  a^exp(-  /k„T)  so  that  a^/c^  =  N  = 

ns  ns  c  B  ns  B  ns  c 

2. 5x10^ (T^/ 300) ^cm  At  high  electric  fields  such  as  those 
encountered  in  the  CCD  fringe  field  region  (  10^V/cm) ,  the  thermal 
emission  rate  increases  only  slightly  which  can  be  largely  accounted 
for  by  lowering  of  the  crap  potential.  The  field  dependence  of 
the  frequency  factor  'a'  is  expected  to  be  weak  [15].  Thus  the  main  hot 
electron  effect  is  the  reduction  of  tie  capture  rate.  This  effect  can 


be  characterized  by  replacing  the  lattice  temperature  by  the  electron 
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temperature,  I.e.  a/c  ■  2.5x10  (T  /T  )  '  .  This  dependence  can 

ns  6  L 

be  understood  in  terms  of  Lax's  cascade  theory  [16].  Using  the 
Maxwellian  distribution,  this  formula  fits  the  experimental  observation 
in  bulk  silicon  quite  well. 

The  electric  field  and  the  carrier  temperature  are  generally 
complicated  functions  of  space  and  time.  The  time  dependence  of  the 
electric  field  near  the  end  of  the  gate  where  the  charge  is  transferred 
to  the  next  gate  is  shown  in  Figure  3.6  for  three  carrier  densities. 


time  [nsi 

Figure  3.6  Electric  field  versus  time  during  transfer 
at  the  end  of  the  transfer  gate.  Gate 
length  is  L  »  2  um 

These  results  show  that  the  electric  field  is  nearly  independent  of  time 

only  for  very  small  carrier  density,  n.  For  this  case,  an  average 

capture  rate  can  be  computed  from 
-1  'L 

c.  »  L  [c  (E)/c  (0)  ]dx  (3.12) 

q  ns  ns 

which  is  plotted  in  Figure  3.7  as  a  function  of  gate  length.  For  large  L, 
c  is  close  to  unity  since  the  electron  temperature  is  high  only  in  a  small 
region  near  the  gap  between  the  gates.  For  small  L,  c  is  less  than  unity. 
The  decrease  in  the  electron  capture  rate  at  small  L  causes  a  reduction 
in  both  the  interface  trapping  noise  and  transfer  inefficiency  due  to 
carrier  trapping  by  interface  states. 
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Figure  3.7  Normalised  average 
caputre  coefficient  versus  the 
gate  length.  Parameters  used 
are  given  in  the  text. 


3 . A  Summary  of  Results 

In  summary,  the  heating  of  the  charge  carriers  in  CCD's  affects 
the  device  performance  as  follows: 

(1)  The  transfer  speed  is  considerably  reduced  because  of  the 
reduction  of  the  mobility  as  already  pointed  out  by  Elsaid  et  al  [A] 
and  shown  in  Figure  3. A.  Diffusion  transfer  becomes  important  only 
near  the  end  of  the  transfer  period. 

(2)  The  noise  and  the  transfer  inefficiency  due  to  carrier 
trapping  by  interface  states  are  reduced  due  to  decrease  in  electron 
capture  rate  which  is  only  important  at  short  gate  length. 

One  of  the  major  assumptions  in  our  calculations  was  that  the 

charge  carrier  heating  is  stationary,  i.e.  the  time  involved  is  longer 

than  the  energy  relaxation  time  which  was  found  to  be  in  the  order  of 
-12 

10  s  both  experimentally  and  theoretically  for  silicon  at  high 
electric  field  [11].  In  Figure  3.3,  we  see  that  for  gates  of  1  „m 
this  condition  is  barely  valid.  The  calculation  for  non-stationary 
heating  is  much  more  involved.  One  expects  a  large  range  of  ohmic 
mobility  and  therefore  faster  transfer  [17]. 


IV. 


TRAPPING  NOISE  IN  CHARGE  COUPLED  DEVICES 


4. 1  Problem  Identification 

Trapping  noise  due  Co  random  transitions  of  electrons  and  holes 
at  interface  states  has  been  investigated  by  several  authors  [18-21]. 

The  mean  square  fluctuation  of  charge  per  unit  area  was  given  by 

•  ’V"TTi“'  (4'1) 

^TTS  t^ie  area^  density  of  interface  states  per  eV  averaged  over  the 

energy  gap.  A  is  the  active  device  area,  q  is  the  elctron  charge 
-19 

(1.6x10  Coulomb),  k  is  the  Boltzmann  constant.  T  is  the  absolute 

B 

temperature.  Equation  (4.1)  is  the  noise  accquired  by  a  charge  packet 

in  a  single  transfer  and  it  is  based  on  the  following  assumptions. 

(i)  Only  thermal  emission  of  crapped  electrons  and  holes  are 
included.  Random  thermal  capture  events  are  ignored. 

(li)  Time  and  space  variation  of  electron  density  during  the 
charge  transfer  are  neglected. 

(iii)  Complete  charge  transfer  is  assumed. 

(iv)  Interface  states  are  assumed  filled  to  the  equilibrium 
level. 

The  last  assumption  was  removed  by  Carnes  and  Kosonocky  [21] 
whose  numerical  solutions  for  high  clock  frequencies  can  be  approximated 


by  the  following  equation. 


<qss>/A  *  qVNITSt“2<q»s”/fcl)  (i-2) 

This  result  has  an  additional  factor,  (c  n/ f  , ) ,  over  the  simple  solution 

ns  cl 

given  by  Equation  (4.1).  This  is  the  fraction  of  the  interface  states 


which  are  filled  when  thermal  emission  begins. 

In  this  section,  the  other  three  neglected  factors  listed  above 


.2 


Theoretical  Analysis 


When  the  assumptions  are  all  removed,  the  simplified  mean 
square  fluctuation  of  charge  per  unit  area,  given  by  Equation  (4.1), 


is  generalized  to  the  following  approximate  form. 

<q^  >/A  =*  qkt,TN__c,F 
'ss  M  B  ITS 

where  the  factor,  F,  has  three  parts  given  by 
F  -  Fi  +  F2  +  F3 
They  are  defined  below. 
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2  Cf  Cf 

d£  f  ensnTSdt’ ^exp (_2/  ensdt ' ) 1 

£  j  C0f  1 


(4.7) 


Here,  £=«E/k  T.  £  is  roughly  the  intrinsic  Fermi  level  below  which  the 

D  2 

interface  state  levels  are  assumed  to  be  all  filled  at  all  time,  c  and 

ns 

e  are  the  thermal  capture  and  emission  rates  of  electrons  at  the 
ns 

interface  state  whose  energy  level  is  E  or  £  =>E  /k„T.  n=n(x,t)  is  the 

s  s  s  B 

instantaneous  electron  concentration  at  position  x.  The  integration  time 
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llmics  are  defined  as  follows  t^  is  the  time  when  the  charge  begins 

to  move  into  the  receiving  ga  e.  t  is  the  time  when  the  charge  under 

on 

the  receiving  gate  exceeds  th<  background  charge  csuch  as  tfac  from  a 

..  .  -1  . 

fat  zero)  already  under  the  receiving  gate  before  t ^ .  cc+  'rcl  ls 

the  time  when  the  charge  unde  *  the  receiving  gate  has  reached  a  maximum 

and  begun  to  decrease  or  spreid  out  to  the  next  gate.  f  ^  is  the  clock. 

frequency.  t  ^  is  the  time  when  the  most  of  the  charges  have  been 

transferred  to  the  next  gate  i rom  the  receiving  gate  such  that  the  charge 

under  the  receiving  gate  has  cropped  to  the  background  level.  t  +(f  ^)  ^ 

is  the  end  of  one  charge  tranc.fer  cycle,  t ^  is  approximately  tc+(f  ^)  ^ • 

The  energy  E  or  Z  =E  ,  k  T  seoarates  the  faster  interface  states 

from  the  slower  interface  stales.  Thus,  the  states  between  the  conduction 

band,  E  =0,  and  E,  have  larger  thermal  emission  rates  e  (E  <E, )  than  those 
C  l  ns  s  1 

between  and  the  intrinsic  i-'ermi  energy,  E0 .  This  energy  boundary 
which  divides  the  shallower  and  faster  states  from  the  slower  and  deep 


states  is  somewhat  arbitrary.  A  rough  definition  is 

e  (E  )  =  c  <n(t  ,) >  (4.8) 

ns  1  ns  or 

where  the  electron  concentration  is  averaged  over  the  gate  at  t=t^-  and 

the  thermal  equilibrium  or  mass  action  law  relationship 

e  (E_)  =  c  (E, )N  exp(-E/k  T)  (-.9) 

ns  1  ns  1  C  B 

is  used.  N  is  the  effective  density  of  states  of  the  conduction  band. 


The  Equation  (4.3)  with  the  three  factors  is  derived  from  the 


rate  equation  of  electron  density  trapped  at  the  interface  state, 

dnTS/dC  "  Cnsn(NTTS  ‘  nTS}  "  enSnTS  +  r3(t)  (i-10) 

where  r^t)  is  a  random  generating  function.  Using  the  small-signal 

expansion  for  small  amplitude  noise,  nTS(x,t)  =  n°g(x,t)  +  in^Cx.t), 
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the  solution  for  the  fluctuating  part  of  the  trapped  electron  density  is 


t’ 


<5nTS(x,t) 


dt ' r  ( t ' )  •  exp 


dt"[c  n(x, t")  +  e  ] 
ns  ns 


(4.11) 


0 


Taking  the  ensemble  average  over  the  active  area  of  the  device,  the 
trapped  charge  fluctuation  is  then  given  by 

eG  Cf 

d€  |  dt[(»TTS  -  n°s)cnsn  +  n^sena]  x 


'0 


x  exp 


-  2 


dt ' (c  n  +  e  ) 
ns  ns 


(4.12) 


where  use  is  made  of 


<rs(tl)rs(C2)>A  "  [(NTTS  -  IlTS) Cnsn  +  Wns^W  (4‘13) 


4.3  Summary  and  Qualitative  Interpretation 

The  third  factor,  F^,  given  by  Equation  (4.7)  reduces  to  the 
simple  result  given  by  Equation  (4.1)  for  constant  and  long 

transfer  time,  i.e. 

F3  =  Ini  (4.14) 

This  is  the  case  where  the  major  contribution  of  noise  comes  from  electron 
trapping  at  the  deeper  interface  levels.  If  the  interface  state  density 
depends  on  energy,  then 

F3  .  l„2  •  tNTTS<Emajt>/NTTS)  (4.15) 

where 

E  -  k  T  [c  N„(t,-t  c)/in2] 
max  B  ns  C  f  of  ‘ 

and  is  the  average  over  the  energy  gap  as  defined  previously  and  it 

cancels  out  the  same  factor  in  Equation  (4.3). 


If  the  electron  density  is  high  or  the  width  of  the  charge  packet  is 

small,  i.e.  if  c  n(t  c-t  )  <  1,  then  Equation  (4.15)  is  diminished 
ns  of  on 

by  this  factor  to  give 


r  =  tn2  •  (N  (E  )/N  ]-c  n(t  -t  ) 

j  TTb  max  TTb  ns  of  on 


(4.16) 


This  differs  from  Equation  (5.2)  slightly  because  electron  capture, 


which  is  incomplete  in  the  time  interval  t  ,-t  ,  can  continue  until 

of  on 

the  charge  packet  has  completely  transferred  out  of  this  gate. 

The  factor  is  associated  with  incomplete  filling  of  the 

interface  state  levels.  When  all  the  levels  are  filled,  n„„(t  )=N^,„_, 

TS  on  TTb 

and  F^=0.  Due  to  incomplete  filling,  the  density  of  the  occupied  and 
unoccupied  interface  states  will  both  fluctuate  over  the  entire  active 


device  area  and  contribute  noise,  but  at  each  interface  state,  the 


densities  are  correlated  since  nTS(x,t)  +  [NTT  (x)  -  nTS(x,t)]  =  NTTS(x). 

The  factor  F^  comes  from  the  filled  traps  and  from  empty  traps. 

The  second  factor,  F  ,  comes  from  the  random  fluctuation  of  the 

capture  transition  into  the  shallow  energy  level  traps  between  E^,  and  E^. 

This  noise  can  become  dominant  when  e  (t  ,-t  )  is  large  and  when  the 

ns  of  on 

electron  density  is  high.  At  low  electron  concentrations,  its  contribution 
may  be  comparable  to  the  noise  due  to  random  emission  and  incomplete 
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V.  SMALL-SIGNAL  IMPEDANCE  OF  INTERFACE  STATES  IN  CCD  OR  MOS  DEVICES 

5 . 1  Problem  Identification 

The  major  cause  of  transfer  inefficiency  in  SCCD  operated  at 

medium  frequencies  is  the  trapping  of  minority  carriers  by  the  interface 

states  located  at  the  Si-Si0o  interface  [22].  The  introduction  of 

background  charge  or  fat  zero  eliminates  most  of  the  "fixed  loss"  due 

to  interface  state  trapping  under  the  gate.  The  major  source  of 

transfer  inefficiency  then  comes  from  interface  states  at  the  edges  of 

the  transfer  electrodes.  This  occurs  because  the  area  of  interface 

over  which  a  charge  packet  is  stored  is  larger  than  that  occupied  by 

background  charge  [23].  This  is  sometimes  referred  to  as  the  non-linear 

1-ss  since  the  loss  is  a  function  of  signal  size  [24].  Another  major 

source  of  transfer  inefficiency  is  the  trapping  of  minority  carriers 

by  the  interface  states  in  the  gap  between  the  electrodes.  This  can 

be  reduced  by  overlapping  the  gates  of  CCD  [25].  The  role  of 

interface  states  in  causing  trapping  noise  was  studied  in  detail  in 

the  previous  sections,  which  showed  that  the  following  interface  state 

parameters  are  required  to  predict  the  noise  and  trapping  inefficiency: 

the  density  of  the  interface  states,  N^^CE),  and  the  thermal  emission 

and  capture  rates  of  electrons  (for  n-channel  CCD),  e  (E)  and  c  (E), 

ns  ns 

all  as  a  function  of  energy.  In  this  section,  we  describe  our  effort 
to  use  the  small-signal  equivalent  circuit  model  to  obtain  these 
interface  states  parameters. 
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5 . 2  Small-signal  Equivalent  Circuit  Models  for  an  MOS  Capacitor  with 
1 nter1 ace  States 

The  Shockley-Read-Hall  (SRH)  model  of  thermal  generation- 

recombinat ion-trapping  at  an  Imperfection  center  has  been  applied  to 

one-level  interface  states  by  Sah  [26].  The  resulting  equivalent  circuit 

is  shown  in  Figure  5.1.  For  the  interface  branches,  Cgg  is  the  trapped 

charge  storage  capacitance.  ^SSp  i-s  c^e  hole  capture  conductance  and 

G  is  the  electron  capture  conductance.  They  are  defined  as  follows: 
ssn  1 


_1_ 

N  „ 

t  (1-f 

kRT 

TTS 

ss 

2 

=  -S__ 

c 

N  N  „„ 

ns 

S  TTS 

2 

=  -9— 
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P  N 

V 

ps 

S  TTS 

ss 


(1-f  ) 

ss 


(5.1) 

(5.2) 

(5.3) 


where  is  the  areal  density  of  a  single-level  interface  state.  f 

TTS  •  ss 

is  the  electron  occupancy  factor  of  these  states.  cng,  c^g  are  the 

3 

electron  and  hole  capture  rate  constants  in  cm  /sec  respectively. 

N  ,  are  the  semiconductor  electron  and  hole  volume  concentration 

in  #/ cm^  at  the  interface  respectively.  Cq  is  the  oxide  capacitance. 

Neglecting  bulk  generation-r ecombination-trapping  centers,  one 

can  simplify  the  circuit  in  Figure  5.1  to  one-lump  equivalence  as 

shown  in  Figure  5.2  for  an  n-type  substrate  device  biased  into 

accumulation  [26] .  C  is  the  electron  storage  capacitance  in  the  bulk. 

n 

It  is  given  by 

r oo 

0  f 


C  =  ,  _ 
n  k3T 


Ndx 


(5.4) 


0 


Similarly,  if  :he  device  is  biased  into  inversion  and  the  interface 
edge  effect  is  neglected,  the  equivalent  circuit  of  Figure  5.1  can  be 


Figure  5.1  The  small-signal  equivalent  circuit  for  an  MOS  capacitor 
deveoped  by  Sah  [26]  for  a  single-level  interface  state 
and  a  single-level  bulk  recombination  center. 


css  Gssn 


Figure  5.2  Simplified  one-lump  equivalent  circuit  for  an 
n-type  MOS  capacitor  with  surface  accumulation 
of  majority  carriers. 
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simplified  to  chat:  shown  in  Figure 
in  the  bulk,  and  is  given  by 


3.3  [26]. 


V  is  diffusion  admittance 
d 


Yd  =  k0T 


D  P  (1  + 
P  B 


1  /■'’ 

:  )  /(D  t  ) 

po  p  po 


T  is  the  minoricv  carrier  lifetime.  ?„  and  D„  are  the  minor i tv 
po  •  BP 

carrier  concentration  and  diffusivity  respectively.  Under  some 

conditions,  Y1  is  much  smaller  than  other  circuit  elements  and  can  be 
d 

treated  as  an  open  circuit. 

In  practice,  interface  states  are  not  single-level  but  are 
distributed  in  energy  md  space.  The  equivalent  circuit  for  interface 
states  distributed  in  energy  can  be  described  by  a  parallel  array  of 


R-C  branches,  one  for  each  energy  as  shown  in  Figure  3.1a  [26,  27,  28]. 

This  array  can  be  simplified  by  integrating  over  bandgap  energies  to  give 

the  series  equivalent  circuit  of  Figure  t.lb.  Asymptotic  formulae  of 

C  and  G  for  <<\  and  wt  >>1  were  derived  in  reference  [28],  where 
se  se  ss  ss 

=  C  / f G  ( 1— f  )]  in  the  case  of  an  n-tvpe  accumulation  region, 

ss  ss  ssn  ss 

At  verv  low  frequencies  with  .  t  •'<!.  we  have 
n  ss 


with  *  r  <l , 
s  s 

c 

se 


we  have 
q  , ^TTS 


(5.6) 


r, 

se 


>5.7) 


Equation  (5.6)  comes  from  the  interface  state  a  i:-  ■.  n  1  v  -nd  is  valid 

for  ail  energy  levels  of  the  interface  states.  1  •  '  re  it  can  be 

used  to  obtain  o  from  experimental  data  as  '.  v:  •  >:  the 

interface  state  branch  can  be  extracted  from  the  data. 

The  effect  of  interface  states  distributed  into  the  oxide  has 
been  studied  using  i  two-step  model,  consisting  of  6RH  transitions  between 


Figure  5.3  Lumped  equivalent  circuit  model  for  an  n-type 
MOS  capacitor  biased  into  inversion  range  but 
neglecting  the  interface  edge  effect  to  be 
illustrated  in  Figure  5.5. 


Figure  5.4  Lumped  equivalent  circuit  model  for  interface 
states  distributed  in  energy.  (a)  Each  series 
branch  corresponds  to  states  at  one  energy. 

(b)  A  one-lump  equivalent  circuit  for  the 
the  distributed  circuit  shown  in  figure  (a). 
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Che  band  and  the  interface  states  and  elastic  tunneling  transitions 

between  the  interface  states  and  the  oxide  traps  distributed  in  the 

oxide.  The  effect  from  this  distributed  tunneling  model  is  most 

significant  on  the  frequency  dispersion  of  MOS  capacitance  and 

conductance  at  low  frequencies,  especially  for  the  equivalent  series 

conductance  G  which  drops  at  low  frequencies, 
se 

The  experimental  data  on  frequency  dispersion  for  the  aevice 
biased  into  inversion  cannot  be  explained  by  the  simple  model 
snown  in  Figure  5.3  since  the  minority  carrier  diffusion  admittance 
is  so  small  that  the  interface  branch  could  hardly  have  any  effect 
on  the  frequency  dispersion  of  the  device  admittance.  Instead,  the 
interface  edge  region  (IER),  as  shown  in  Figure  5.5,  must  be  included 
in  the  model.  A  circular  geometry  is  used  to  simplify  the  analysis. 

A  two-dimensional  equivalent  circuit  model  is  formed  by  adding  a 
transverse  branch  (in  the  radial  direction)  which  accounts  for  the 
int  irface  edge  effect.  This  is  shown  in  Figure  5.6a.  The  z  direction 
represents  the  usual  MOS  capacitor  branch  as  in  Figure  5.3.  The 
transverse  or  r  direction  represents  the  IER  branch  which  is  essentially 
like  the  usual  MOS  capacitor  branch  but  with  radially  distributed  inter¬ 
face  states  in  the  space  charge  region  of  the  IER. 

The  circuit  elements  for  the  transverse  branch  are  defined  in 

Table  5.1.  GC  is  treated  as  a  short  while  C  ,  G  and  GC  are  treated 
p  n  ssn  n 

as  open  because  the  semiconductor  surface  under  the  gate  is  inverted. 

The  resulting  circuit  is  shown  in  Figure  5.6b  which  is  further  simplified 

to  Figure  5.6c  since  G  is  usuallv  much  larger  than  uC  and  since  C 

ssp  ss  p 

in  weak  inversion  is  much  smaller  than  C  .  The  circuit  in  Figure  5.6c 

ss 


Figure  5.5  MOS  capacitor  structure.  (a)  The  cross-sectional 
view  showing  the  interface  edge  region  (IER)  where 
interface  recombination  may  dominate,  (b)  The  top 
view  of  the  interface  edge  region. 
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Table  5.1  Definition  of  circuit  elements  in  the  transverse  branch 
of  Figure  5.6 


ss 

„t 

Jssn 


series  interface  s  ate  capacitance  of  IER 


series  conductance  of  interface  states  in  IER  due  to 
electron  capture-emission 


ssp 


series  conductance  of  interface  states  in  IER  due  to  hole 
capture-emission 


space  charge  layer  capacitance  of  IER 


minority  carrier  diffusion  admittance  along  the  transverse 
direction 


conductance  of  holes  in  valence  band  from  inverted  region 
to  IER  along  the  transverse  direction 


conductance  of  electrons  in  conduction  band  from  inverted 
region  to  IER  along  the  transverse  direction 


is  then  converted  to  Figure  5.6d  where 
,t 


c  =  c  +  c 

s  w  w 
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(5.8) 

(5.9a) 


(5.9b) 


^SP  ^  ^ 

(c)  (d) 

Figure  5.6  Complete  and  simplified  equivalent  circuit  models  for 
n-type  MOS  capacitor  in  inversion  with  interface  edge 
effect.  (a)Two-dimensional  model.  (b),(c),(d)  simplified 
models  where  the  diffusion  admittance  Y,  and  Y1:  are  opened 


For  high  and  low  frequencies,  we  have 
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Note  that  the  correction  C'  (w+0)  =  C  +  C  must  be  made  for  a 

ss  ss  p 

strongly  inverted  surface  when  is  comparable  to  or  greater  than 
C  .  The  value  of  the  circuit  elements  can  be  extracted  from  the 

5  S 

curves  of  calculated  Cgg  and  Ggg/iu  versus  frequency  from  experimental 
data,  using  the  asymptotic  formulae.  Finally,  we  note  that  the  interface 
edge  effect  is  not  important  Ln  the  accumulation  region  since  the 
majority  carrier  conductance  ind  storage  capacitance  is  so  large  as 
to  short  out  the  IER. 


5.3  De vice  Fabrications  and  Measurements 

The  devices  used  in  this  study  are  made  on  n-tvpe  silicon  with 
<111>  surface  orientation.  The  geometry  of  the  devices  was  shown  in 
Figure  5.5.  The  diameter  of  the  capacitor  is  30  mils.  Phosphorus 
predeposition  was  used  to  getter  impurities  and  to  form  an  n+  layer  for 
good  back  contact.  The  oxide  was  grown  thermally  at  1200°C  in  dry 

oxygen,  followed  by  20  minute  annealing  in  dry  argon  ambient.  The 

o 

thickness  of  oxide  is  about  2600  A.  Room  temperature  bias  stress  tests 
were  performed  and  the  1  MHz  C-V  curves  showed  negligible  drift  of  ions 
m  the  oxide.  This  is  important  since  the  shift  in  C-V  curves  due  to 
mobile  ions  will  cause  significant  error  in  the  total  density  of  inter¬ 
face  states  as  well  as  in  the  frequency  dispersion  data. 
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The  actual  substrate  doping  densities  of  the  device  was  obtained 

from  low  temperature  equilibrium  C-V  curves  with  a  fast  negative 

-2 

voltage  ramp.  C  was  then  plotted  against  V  and  least-square  fitted 

(j 

-2 

to  obtain  N  _  =*  2/qe  (dC  /dV.) .  The  minority  carrier  lifetime,  t  , 

DD  s  G  po 

was  measured  to  be  7.4  usee,  using  the  capacitance  transient  method  of 

Sah  and  Fu  [29].  The  minority  carrier  diffusion  admittance  in  Equation 

(5.5)  is  Y,  =  1.73x10  ^  +  joj  6.43x10  ^  mho/cm^  when  cut  <<  1.  Since 
d  po 

Y,  is  small,  the  approximation  that  Y,  could  be  treated  as  an  open 
d  d 

circuit  in  the  equivalent  circuit  compared  with  other  elements  is 
justified. 

Three-terminal  capacitance  and  conductance  measurements  were 

made  from  1  Hz  to  10'>  Hz  using  a  PAR-186  phase-sensitive  detector- 

amplifier  connected  to  a  Hewlett-Packard  Multiprogrammer  .utomatic 

data  acquisition  system  and  a  Hewlett-Packard  1000  minicomputer  system. 

The  block  diagram  of  the  measurement  system  is  shown  in  Figure  5.7. 

The  details  of  the  admittance  bridge  are  shown  in  Figure  5.8.  The 

transfer  function  is  proportional  to  G  +  jwC  if  (G  +  jwC  )  <<  R  . 

mJm  m  m 

The  automated  data  acquisition  system  is  controlled  by  a  user  program 
which  steps  the  bias  voltage  applied  to  the  device  and  waits  long 
enough  for  output  of  the  lock-in  amplifier  to  settle.  It  then  takes 
a  predetermined  number  of  readings.  The  data  are  averaged  and  stored 
in  a  magnetic  disk  file  for  further  analysis.  Very  accurate  experimental 
data  are  obtained  this  way. 

Accurate  phasing  for  conductance  is  achieved  by  switching-in  a 
lossless  capacitor  with  a  value  of  the  order  of  the  oxide  capacitance  and 
then  looking  for  the  phase  at  which  there  is  no  change  in  conductance  readin 
when  the  lossless  capacitor  is  switched  in  and  out.  Accurate  phasing 


Figure  5.8  The  circuit  diagram  of  the  admittance  bridge  using 
to  measure  the  MOS  admittance  from  1  Hz  to  100  KHz 
using  the  PAR-186  phase-sensitive  detector,  e^  is 
the  reference  audio  signal  and  eQ  is  proportional  to 
a.c.  current  flowing  in  the  MOS  admittance  which  is 
connected  to  the  differential  input  of  the  PAR-186. 

Cs  is  the  standard  calibration  capacitance  and  Cm  and 
Gm  are  the  equivalent  parallel  capacitance  and  conductance 
of  the  MOS  capacitor  under  test. 
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for  conductance  measurement  is  achieved  in  a  similar  manner  by  switching- 
in  a  deposited  precision  carbon  resistor  with  very  low  capacitance. 
Calibration  of  the  capacitance  measurement  is  made  by  using  a  Boonton 
71-3A  push-button  precision  capacitor  which  has  been  calibrated  at  each 
measurement  frequency  by  capacitance  bridges.  Calibration  of  the 
conductance  measurement  is  done  by  taking  advantage  of  the  fact  that 
G/uj  is  just  an  equivalent  capacitance.  The  stray  capacitance  and 
conductance  of  experiment  jigs,  cables  and  header  were  subtracted  out 
from  the  readings  of  the  data  acquisition  program. 

1  MHz  capacitance  measurement  was  made  by  a  computer  controlled 
daca  acquisition  system  using  a  Boonton  72A  capacitance  meter.  The 
amplitude  of  the  small  ac  signal  used  in  this  study  was  between  1  to  10 
mV  whic  i  is  small  compared  with  the  thermal  voltage  k  T/q  =  26mV,  so  that 

D 

the  linear  small-signal  model  is  valid  and  the  assumption  that  a  constant 

density  of  states  is  being  probed  hv  small  ac  signals  is  justified. 

The  experimental  frequencv  dispersion  curves  of  capacitance  and 

conductance  versus  bias  for  a  typical  device  are  shown  in  Figures  5.9 

and  5.10  respectively.  Appreciable  frequency  dispersion  was  observed 

in  the  frequency  range  from  1  Hz  to  10^  Hz.  Series  equivalent 

circuit  elements  C'  ,  C'  rel  ited  to  interface  states  are  calculated 
ss  ss 

from  the  measurements  of  C  and  G  bv  transforming  the  circuit  in 

mm-  s 

Figure  5.11(a)  to  the  circuit  in  Figure  5.11(b).  In  doing  so,  we  have 

assumed  that  the  surface  potential  at  a  given  bias  is  the  same  for 

different  frequencies.  Thus,  the  semiconductor  capacitance  C  can  be 

s 

obtained  from  the  1  MHz  C-V  curve.  C  ,  the  oxide  capacitance  was 
also  calculated  from  the  1  MHz  C-V  curve  in  strong  accumulation 
using  the  method  of  McNutt  and  Sah  130],  Accurate  determination 
"f  Che  oxide  capacitance  is  important  in  this  studv 


CAPACITANCE  C  (pF) 


D.C .  APPLIED  GATE  VOLTAGE  (V) 


Figure  5.9  The  measured  equivalent  parallel  capacitance  of 
MOS  capacitor  24-12C-H6  as  a  function  of  d.c. 
applied  gate  voltage  with  signal  frequency  as 
the  parameter. 


CONDUCTANCE 


Figure  5.10  The  measured  equivalent  parallel  conductance  of 
MOS  capacitor  24-12C-H6  as  a  function  of  d.c. 
applied  ;ate  voltage  with  signal  frequency  as 
the  paraneter. 


Figure  5.11  Equivalent  circuit  model  of  a  MOS  capacitor. 

(a)  The  measured  equivalent  parallel  admittance 
and  (b)  the  equivalent  circuit  applied  to  both 
surface  accumulation  and  inversion  conditions. 


because  in  Che  accumulation  range,  the  interface  states  capacitances 

and  the  majority  carrier  storage  capacitances  .ire  so  large 

that  the  measured  canacitance  C  is  close  to  the  oxide  capacitance. 

m 

In  accumulation,  we  have  C'  =  C  and  G'  =  G  as  in  Figure  5.  A 

ss  se  ss  se 

and  5.2,  In  inversion,  C'  and  G'  are  given  bv  Equation  (5.9)  and 

ss  ss 

(5.10). 

5 . A  Discussion  of  Results 

C'  and  C'  /w  versus  frequency  are  plotted  in  Figure  5.12  and 
ss  ss 

5.13  for  Device  2A-12C-H6  biased  into  accumulation  and  inversion 

respectively.  From  Figure  5.L2  and  5.13,  we  note  that  Cr.e  tsymptoti- 

formulae  G'  (w-  ”)  -*•  G  C  and  C'  ( ■*  =  0)  -*•  C  indeed  noli;.  ~'.e  detisif 
ss  ssp  ss  ss 

of  interface  states  at  this  potential  can  be  calculated  using 

0  ? 

Equation  (5.6),  P„  -  C  / q“  -  C'  / q“,  which  applies  to  both  accumulation 
and  inversion  since  the  equivalent  circuits  have  the  same  form 
(Figures  5.2  and  5.6d).  The  density  of  interface  states  was  also 
obtained  from  Tertian's  method  [31],  l.e.  by  differentiating  the  1  MHz  C-V 
curve,  dC/dV.  The  two  results  agree  with  each  other  fairly  well  in  the 
inversion  region  (V  -  V  <  0)  as  indicated  in  Figure  5.1A.  In  the 

S  r  ^ 

accumulation  region,  the  7  ^  from  low  frequency  capacitance  is  2  to  3 

times  greater  than  that  obtained  by  Taman' s  technique. 

The  experimental  accuracy  of  the  capacitance  measurements  in 

A 

inversion  for  frequencies  higher  than  10  Hz  is  not  good.  This  fact  is 
clearly  shown  in  Figure  5.9.  The  capacitance  in  the  inversion  region 
increases  only  slightly  for  frequencies  above  1C)14.  The  range  of  surface 
potentials  available  in  this  study  is  limited  bv  experimental  accuracy 
as  well  as  by  the  pinning  of  surface  potential  due  to  high  density  of 
interface  states. 


1  10  102  103  104  105 


FREQUENCY  (Hz) 

Figure  5.1 1  and  G^g/'a  versus  frequency  of  MOS  capacitor  24-12C-H6 

at  two  interface  state  energies:  E. +127.1  and  E.+130.3  mV. 
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Figure  5. 13  C'  and  O'  /oj  versus  frequency  for  Device  24-12C-H6  at  a 
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interface  state  level  110.6  mV  above  the  intrinsic  Fermi  level 
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The  behavior  of  Che  low  frequency  C-V  in  Figure  5.9  can  be 

readily  explained  using  Che  equivalenc  circuic  of  Figure  5 . 6d . 

5  2  4  2 

C  is  escimaced  Co  be  8.5x10  pF/cm  from  Figure  5.13.  Cg =  10  pF/cm 

4  2 

and  -  1.33x10  pF/cm  .  Therefore  as  ^  -*■  0,  C$s  in  Figure  5 . 6d 
dominaces  and  Che  measured  capaciCance  approaches  Che  oxide  capacicance 
as  shown  in  Figure  5.9. 

5 . 5  Summary  of  ResulCs 

Simple  one-lump  equivalenc  circuiCs  were  derived.  They  were 
found  co  be  quice  successful  in  explaining  Che  frequency  response  of 
incerface  scaces.  '.'Then  Che  device  is  biased  inco  accumulacion,  Che 
series  equivalenc  circuic  for  incerface  scaces  cogecher  wich  a  cwo-scep 
model  involving  Shockley-Read-Hall  CransiCions  and  Cunneling  processes 
works  well.  In  Che  case  of  an  inverced  surface,  a  model  which  Cakes 
inco  accounc  the  incerface  edge  effect  is  proposed  and  found  co  be  in 
good  agreemenC  with  experimental  data.  Both  circuic  models  allow  one 
co  calculate  the  densicy  of  incerface  states  from  low  frequency  inter¬ 
face  state  capacitance  and  Che  results  were  compared  wich  those  obtained 
from  Terman's  high-frequency  C-V  method  [31].  The  agreement  between 
these  results  was  quice  good.  These  equivalenc  circuits  contain  only 
linear  elements  (capacitances  and  conductances).  The  interface  state 
properties  can  be  obtained  from  the  experimental  values  of  the 
:ircuit  elements  from  the  low  and  high  frequency  measurements  of  MOS 
capacitance  and  conductance. 
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VI.  SULK  TRAPS  IS  BURIED-CHANNEL  CHARGE-COUPLED  DEVICES 
o . 1  Problem  Identif ication 

In  buried-channel  charge-coupled  devices  (BCCD),  Che  charge  is 
stored  away  from  Che  Si-SiO?  interface  so  that  carrier  trapping  at 
the  interface  states  can  be  neglected.  In  this  case,  the  trapping  at 
bulk  impurities  and  defects  limits  the  device  performance.  One  may 
introduce  a  small  background  charge  (or  fat  zero)  to  reduce  the  transfer 
inefficiency  for  a  BCCD  in  a  similar  way  to  that  used  for  a  surface 
channel  CCD.  In  the  presence  of  background  charge,  edge  effects  are 
the  dominant  loss  mechanism  because  the  signal  charge  occupies  a 
larger  volume  than  the  background  charge. 

Unlike  the  interface  states,  the  bulk  states  (or  bulk  traps) 
have  discrete  energy  levels  in  the  energy  gap.  The  thermal  emission 
rate  of  electron  trapped  by  bulk  states  is  related  to  the  activation 
energy  and  temperature  by  the  Arrhenius  equation 


e 

n 


i  v  N  esp  (-E  ,  T>  . 
n  t  c  A3 


Here  7  is  the  capture  cross  section.  v  is  the  thermal  velocitv. 
n  t 

N  is  the  density  of  states.  E  is  the  activation  energy.  We  note 

C  A 

that  the  emission  rate,  e  or  emission  time  constant,  ~  =  1/e  ,  is  in 

n  n  n 

exponential  function  of  temperature.  Each  species  of  bulk  traps  has 
a  characteristic  lifetime  given  by  Ecuation  (6.1). 

Experimental  data  of  transfer  inefficiency  versus  clock 
frequencies  reveal  severaL  peaks  [32,  33].  Similar  peaks  were  also 
reported  in  BCCD  noise  measurements  [34].  It  was  shown  that  for  a 
three  ohase  CCD,  the  transfer  inefficiency  due  to  bulk  state  trapping 
reaches  the  maximum  value  when  the  clock  frequency  is  related  to 
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emission  lifetimes  by  [35] 

f  =  (N  -r  1/3)  [’  ln(3N 
max  z  e  z 


2)1- 


(6.2) 


Here  Is  the  number  of  zeros  preceding  the  first  signal  "one." 

10  12  -3 

Bulk  trap  densities  between  10  and  10  cm  have  been  reported 
[32-34],  Some  of  these  bulk  traps  come  from  metallic  impurities  such 
as  gold,  sulfur  and  iron  unintentionally  introduced  during  processing. 
These  impurities  were  identified  by  matching  the  measured  activation 
energy  with  the  published  data  of  known  bulk  traps  [32].  A  more 
reliable  identification  is  to  match  the  measured  emission  tine  constants 
from  the  BCCD  in  a  wide  range  of  temperatures  with  published  data  of  a 
known  bulk  trap  such  as  gold  [36]  and  sulfur  [37].  The  fact  that  the 
bulk  traps  observed  in  BCCD  could  also  come  from  orocess-induced  defects 
[33,  39]  has  not  been  noted  in  the  literature. 

Ion  implantation  is  ised  to  fabric  r.e  BCCD .  A  large  numoer 
bulk  trap  levels  were  observed  in  silicon  after  ion  implantation  [40]. 

It  is  likely  that  some  of  these  bulk  states  may  still  be  present  after 

9 

high  temperature  annealing  since  bulk  trap  densities  in  the  range  of  10 
to  10^  cm  ^  is  sufficient  to  cause  problems  in  3CCD.  These  bulk  traps 
can  also  come  from  thermally-: ndueed  defects  or  quenched-in  centers. 

Yau  and  Sah  [38]  reported  two  trapping  levels  located  at  0 .  2o4  and 
0.542  eV  below  the  conduction  band  for  n-type  silicon.  Collet  [32] 
observed  two  levels  at  0.25  and  0.54  eV  below  the  conduction  band  in  his 
n-channel  3CCD.  These  could  very  well  be  the  same  quench- in  levels 
reported  by  Yau  and  S.ah. 

The  objective  of  this  project  is  to  accurate 1 v  haracter  i.:o  '  '.e 
buLx  traps  which  may  be  introd  iced  luring  BCCD  processing  steps  such  i- 
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ion-implantation,  expitaxial  growth,  thermal  quenching,  chemical 
cleaning  and  diffusion. 

6 . 2  Experimental  Techniques 

Transient  capacitance  techniques  [41,  42]  are  used  because  of 
their  high  sensitivity  and  also  because  the  experimental  conditions 
of  diode  fabrication  and  measurements  are  very  similar  to  the 
operating  conditions  of  the  BCCD. 

Our  measurement  has  been  computer i?ed .  A  block  diagram  of  the 
apparatus  used  is  shown  in  Figure  6.1.  The  device  to  be  measured  and 
the  temperature  controller  are  shown  in  Figure  6.2.  The  detail 
operation  of  the  temperature  controller  was  described  in  a  previous 
report  [43].  The  temperature  controller  is  implemented  in  software 
using  proportional  plus  integral  control.  The  controller  is  capable 
of  controlling  the  sample  temperature  to  within  +0.1  K  during  the 
measurement  in  the  temperature  range  from  85  K  to  300  K.  The  controller 
controls  the  heater  power  by  altering  the  output  voltage  of  the 
HP  6224B  power  supply  through  the  D/A  card.  The  device  temperature  is 
monitored  by  the  Data  Precision  5-1/2  digit  voltmeter  through  two 
isolated  input  cards. 

The  data  sampling  is  controlled  by  a  separate  program.  The 
sampling  program  starts  an  experiment  by  sending  a  pre-determined  bit 
pattern  to  the  TTL  output  card  which  Ln  term  triggers  the  pulse 
generator  (HP  8010A)  to  generate  a  pulse  of  preset  amplitude  and 
duration.  The  resulting  capacitance  transient  due  to  the  pulse  is 
measured  and  converted  into  digital  torm  by  a  second  5-1/2  digit  DVM. 

The  time  at  which  each  capacitance  data  is  tiken  is  obtaine  i  from  the 
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Figure  6.1  Block  diagram  of  Che  computer  controlled  transient 
capacitance  measurement  setup. 


junction 
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real-time  clock  in  the  minicomputer.  The  capacitance,  time  and  the 
temperature  of  each  data  point  are  stored  in  a  disk  file  for  plotting 
and  further  analysis.  A  typical  capacitance  versus  time  plot  is  shown 
in  Figure  6-3.  The  temperature  during  the  experiment  is  shown  on  the 
upper  part  of  the  plot.  The  small  random  variation  in  temperature  comes 
from  electrical  noise  picked  up  by  the  thermocouple  wire  which  shows 
that  the  device  temperature  was  kept  constant  to  better  than  +0.1  K. 


6.3  Measurement  of  Bulk  Trap  Parameters 


Thermal  emission  rate  is  obtained  by  least-square  fitting  the 


C-t  (capacitance  versus  time)  data  in  Figure  6.3  to 


C(t)  -  C  -  AC-e-C/T*.  (6.3) 

The  thermal  emission  rates  or  the  emission  time  constants 
determined  at  many  temperatures  are  then  fitted  to  the  equation 


t_1  =  eC  =  A  (T/300)m  exp  (-E,/k_T).  (6.4) 

e  n  .\  3 

Equation  (6.4)  is  identical  to  Equation  (6.1)  except  that  the  temperatur 
dependence  of  the  pre-exponential  factor  is  explicitly  displayed 
in  (6.4).  In  general,  the  fitted  thermal  activation  energy  or  the 
energy  level  of  the  bulk  traps  depends  on  the  choice  of  m  in  Equation 
(6.4).  The  total  capacitance  change,  AC,  during  the  entire  transient 
is  related  to  the  ratio  of  the  bulk  trap  density  to  majority  impurity 
dopant  or  majority  carrier  concentration  as  follows  [42] 


ic/c  =  m 


'  V 


Equation  (6.5)  is  valid  if  ?.C<'C  or  N_  <<(N_ _  -  M,,).  Thus,  the 

”  XT  DU  AA 

ipproximate  density  of  bulk  traps  can  be  obtained  from  C-t  data  and 
Equa  t ion  (6.5). 
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It  is  usual Lv  desirable  Co  make  a  quick  survey  Co  determine 

it  there  are  bulk  traps  present  in  i  given  sample.  This  Surrey,  which 

we  denote  VSCTS  '.Voltage  Stimulated  Capacitance  Transient  Spectroscopy), 

can  be  performed  using  the  compucer  controlled  data  acquisition  system 

shown  in  Figure  6.1.  In  the  VSCTS  measurements,  the  device  temperature 

is  raised  from  77K  to  300K  either  by  a  programmed  slow  linear  or 

nonlinear  temperature  ramp  (known  as  the  DLTS  method  [44])  or  by 

programmed  small  temperature  steps  (VSCTS) .  Instead  of  sampling  the 

entire  capacitance  transient  at  each  constant  temperature,  the 

capacitance  is  sampled  at  two  preselected  and  fixed  times,  t^  and  tr 

This  is  illustrated  in  the  inset  of  Figure  6.4.  The  capacitance 

difference  at  these  two  times,  C?-C^,  is  then  stored  and  plotted  as  a 

function  of  temperature.  Its  magnitude  is  a  function  of  temperature  or 

the  emission  time  constant,  r and  the  sampling  times  t,  and  t , .  A 

typical  example  from  a  gold-doped  silicon  n+/p  diode  is  shown  in 

Figure  6.-.  The  occurence  of  the  capacitance-d  if  f  erenc.e-sigr.al  peak 

at  a  particular  temperature  (in  this  case  140K  for  the  large  peak)  can 

be  understood  from  Figure  6.4.  Consider  a  fixed  pair  of  sampling 

times,  t^  and  t,,  for  all  temperatures  between  77K  and  300K.  At 

low  temperatures  when  the  thermal  emission  time  constant,  x  ,  is  large 

compared  with  the  sampling  time  interval,  t^-t^,  the  capacitance 

transient  decay  is  slow  and  C,-C  is  small.  As  temperature  rises  and 

x  decreases,  the  caoacitance  transient  decav  rate  increases  and 
e 

C,-C,  will  increase.  However,  at  still  higher  temperatures, 
becomes  so  short  that  the  transient  is  nearly  over  at  t,,  so  that 
C,-C  would  again  be  small.  If  the  capacitance  transient  is  truly 
exponential,  it  can  be  easily  shown  taat  the  thermal  emission  rate  is 
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given  by  [44] 

e C  =  ( t ,  -  c  )  1  •  in  ( t .  / 1 ^ )  (6.6) 

figure  6.4  shows  two  p.  aks  at  140K  and  230K  due  to  the  two 
gold  levels  in  silicon.  The  40K  is  due  to  hole  emission  from  the 
lower  donor  level,  E^+355  mV.  The  higher  cemperature  and  smaller  peak 
at  230K  is  due  to  electron  emission,  even  though  the  silicon  substrate 
is  p-type.  This  is  due  to  the  fact  that  the  emission  time  constant 
depends  on  the  emission  rate  of  both  electrons  and  holes  trapped  at 
a  given  level  [42] 


T  =  l/eC  =  l/(e  +  e  )  =  1/e  (6.7) 

e  n  p  n 

since  at  this  gold  acceptor  level  near  the  midgap  of  silicon,  the  electron 

emission  rate  is  much  higher  than  hole  emission  rate,  i.e.  e  >>e  . 

n  p 

The  smaller  amplitude  of  the  230K  peak  compared  with  the  140K 
peak  is  due  to  the  fact  that  most  of  this  midgap  level  is  occupied  by 
holes  in  a  p-type  silicon  and  the  trapped  electron  density  is  very  small 
comparer  with  the  gold  density.  It  is  approximately  given  by 


U_  =  N’  e  /  (e  +e  )  '  N,  (e  /e  )<<N 
T  TT  p  n  p  Au  p  n  Au 


The  larger  amplitude  of  the  140K  peak  corresponds  to  the  trapped  hole 

concentration  at  the  lower  donor  level  and  it  is  almost  equal  to  the 

gold  concentration,  or  N,  . 

TT  Au 

The  difference  in  capacitance  signal  at  the  two  levels 
provides  another  means  of  determining  the  minority  carrier  as  well  as 
the  majority  carrier  emission  rates  in  one  experimental  transient 


capacitance  run. 
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Figure  6.3  A  typical  computer  output  of  the  capacitance  transient. 

The  upper  figure  is  the  device  temperature  and  the  lower 
figure  is  the  capacitance  transient. 
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In  interpreting  VSCTS  data,  it  is  n.-cessary  to  eliminate 


os 


b .  -* 


signal  from  trapping  in  a  surface  channel.  One  example  of  this  is 
shown  in  Figure  o.3.  This  data  is  taken  from  the  same  device  from 
which  data  in  Figure  6.4  is  taken  except  that  a  MOS  guard  ring  around 
the  n+  contact  is  floated  in  Figure  6.5  while  Figure  6.4  is  obtained 
with  a  negative  voltage  applied  to  the  MOS  guard  ring  to  pinch  off  the 
surface  channel.  The  extra  peak  in  Figure  6.3  comes  from  the  surface 
channel  which  is  pinched  off  in  Figure  6.4. 

Another  example  of  a  surface  effect  is  observed  in  an  aluminum 
on  n-type  silicon  Schottky  barrier  diode.  The  VSCTS  data  taken 
shortly  after  fabrication  is  shown  in  Figure  6.6.  There  are  two  peaks 
at  120  and  300K.  as  well  as  large  noise.  The  300K.  (midgap  leveL)  peak 
and  the  noise  are  related  to  the  high  reverse  leakage  current  of 
the  device  which  is  6 uA  at  -10  volts  at  room  temperature.  Since  the 
n-type  surface  area  surrounding  the  30  mil  diameter  circular  diode  is 
not  passivated,  it  is  likely  that  the  large  leakage  current  is  caused 
by  water  vapor  around  the  edge  of  the  diode  [45’,.  This  is  confirmed 
by  baking  the  diode  at  200C  for  15  minutes  in  dry  nitrogen  ambient. 

The  room  temperature  leakage  current  after  bakeout  was  reduced  to  less 
than  0.5  .jA  at  -10  volts.  The  VSCTS  was  then  taken  and  is  shown  in 
Figure  6.7  where  the  30CK  peak  is  eliminated  and  the  noise  is 
reduced  to  0.1  mpF  which  is  the  instrument  noise  from  the  Boonton  72A 
one  MHz  capacitance  meter. 

n  .  3  Fee  or.  binati  jr  j  Centers  l:  t _ I  an  Implanted  Silicon 

Some  ore!,  iminarv  results  is  ir.g  the  •••nputer  controlled 
an  p;  itunce  transient  measurement  .setup  have  been  obtained  :  r  -  o 
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oxygen  implanted  n+/p  diode  (implanted  at  AFCRL)  previously 

fabricated  in  our  laboratory  [40] .  The  implant  was  done  at  high 

13  2 

energy  (2MeV)  and  high  dose  (10  ions/cm  ).  The  VSCTS  data  for  three 
different  filling  pulse  amplitudes  are  shown  in  Figure  6.8.  The  d.c. 
bias  voltage  was  10V.  Seven  peaks  are  identified.  They  are  labelled  PI 
tnrougn  ?/.  The  relative  change  in  the  size  of  peaks  at  different 
pulse  amplitudes  indicates  that  the  bulk  traps  species  have  quite 
different  spatial  distributions.  The  spatial  distribution  of  each 
bulk  trap  species  can  be  measured  using  the  differential  bias  voltage 
technique  [46],  The  VSCTS  data  of  another  ion  implanted  diode  from 
the  same  wafer  but  annealed  at  400C  for  30  seconds  in  forming  gas 
ambient  is  shown  in  Figure  6.9.  Peaks  P2,  P3  and  P6  are  eliminated; 

P5  is  barely  detectable.  PI  decreases  in  size  while  P4  and  P7 
increase  quite  substantially.  This  result  shows  a  reverse  anneal 
effect  in  which  P4  (160K)  and  P7  (290K)  have  increased  their 
concentration  during  heating  at  4Q0C. 


AL  ON  N-SILICON  SCH0TTKT  BARRIER 
AFTER  BAKEOUT  AT  200  C  FOR  15  MINUTES 
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Figure  6.7  VSCTS  of  the  Al/n-Si  Schottkv  barrier  diode  shown  in 

Fig. 6. 6  after  bakeout  at  200C  for  15  min.  in  dry  nitrogen. 
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;ure  6.3  Voltage  Stimulated  Capacitance  Transient  Spectra  of  an  oxvgen 
(2  MeV)  implanted  n+/ p  Si  diode  after  heated  at  70C  for  30  sec 
during  soldering  of  the  device  dice  to  an  8-pin  7C-3  header. 
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Figure  6.9  Voltage  Stimulated  Capacitance  Transient  Spectra  of  the  oxvgen 
implanted  n+/p  Si  diode  shown  in  Fig.b.S  after  annealing  at 
4ooC  for  30  secs  in  forming  gas. 
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VII.  PERFORMANCE  LIMITED  BY  INTERBAND  IMPACT  IONZATION 

The  longitudinal  electric  field  (longitudinal  is  the  direction 

along  which  the  charges  are  transferred,  i.e.  the  direction  parallel  to 

the  oxiie-sil icon  interface)  in  a  CCD  varies  from  a  minimum  at  the 

cencer  of  a  gate  to  a  maximum  at  the  center  of  a  gap.  For  this 

discussion,  the  magnitude  of  these  electric  fields  are  denoted  by 

E  .  and  E,, 
min  Max. 

In  order  to  achieve  the  highest  possible  charge  transfer  speed, 
E  .  should  be  made  as  large  as  possible  since  the  charge  trasfer  time 

mm  a  r  a 

under  the  gate  is  approximately  giver,  by  the  gate  length  divided  by  the 


drift  velocitv: 


L/v  »  L/tE  , 
d  min 


where  c=u(E  ,  )  is  the  field  dependent  drift  mobilitv.  A  similar  estimate 
min 

can  be  made  for  the  charge  transfer  time  across  the  gap, 

t  =>  d/ v  =  d/ uE  (7.2) 

gap  d  Max 

where  d  is  the  length  of  the  gap  and  =  j(E  ).  If  the  gate  length  is 

cltiX 

large,  the  transfer  time  under  the  gate  may  become  diffusion  limited 

instead  of  drift  limited  as  given  by  Equation  (7.1).  The  gate  transit 

time  for  the  diffusion  limited  condition  is 

t'  =  L2/D  =  (L/uE) (ELq/k  T)  i  (L/uE)*F  (7.3) 

gate  3  gate 

which  is  larger  than  the  drift  limited  case,  Equation  (7.1),  by  a  factor 

EL/(kBT/U  3  Fg,te  <7'4) 

before  velocity  saturation  sets  in  at  high  electric  field. 

It  is  evident  that  the  minimum  electric  field  under  the  gate 


should  be  made  as  high  as  possible  to  obtain  high  transfer  speed. 
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E  .  =  5(z  V/L~) 

min  ox 


has  been  estimated  and  was  given  by  [47] 

(7.5) 

where  V  is  the  d.c.  voltage  difference  between  two  adjacent  gates. 

For  a  CC1)  with  overlapping  gate  or  zero  gap  width,  the  maximum  gap 
field  is  approximately  given  by 


EL,  V/5z 
Max  ov 


(7.6) 


which  must  be  smaller  than  the  breakdown  field  of  silicon,  E  .  Thus, 

DU 


the  minimum  field  must  be  less  than 


E  .  <  25  (z  /L)  *"E 

min  ox  Bb 


(7.7) 


The  gate  transfer  time  from  Equation  (7.1)  using  Equation  (7.7)  is  then 

L  >  L 


t  _  >  L3/(25z2  uF„ _ ) 

gate  ox  uS 


sat 


L/v 


L 


(7.8) 

(7.9) 


sat  sat 

Here,  L  is  the  gate  length  below  which  drift  velocity  saturation 
sat 

occurs  due  to  hot  electron  effect  while  the  gate  voltage  difference,  V, 
is  kept  as  high  as  possible  but  just  below  the  silicon  breakdown  field, 

E__.  Thus,  using  v  =uE  .  ,  L  is  then 
BS  sat  mm  sat 


L  „  =  5z  >jE  _/v 
sat  ox  BS  sat 


(7.10) 

The  adjacent  gate  voltage  difference  required  to  get  the  highest  E 
under  the  gate  before  Si  bulk  breakdown  occurs  from  Equations  (7. 5^  and 
(7.7)  is 

V  5z  E  (7.11) 

ox  BS 

These  simple  physical  considerations  give  us  the  necessary 
design  equations  to  determine  the  ultimate  transfer  time,  which  is 
limited  by  the  transfer  time  under  the  gate.  As  a  numerical  example, 
let  us  assume  that  E  =103  V/cm  above  which  excessive  interband  electron- 

DU 


hoLe  pair  generation  would  destroy  the  zeros.  Let  us  take  z  =0.1  ..m. 


an  electron  surface  mobility  of  p=700  cm“/Vs  and  a  saturation  velocity 

of  v  =4x10^  V/cm.  Then, 
sat 

L  =*  2.09  pm  (7.12) 

sat 

For  wide  gate,  such  as  L=10  pm,  we  have  the  following  results. 

E  .  «  250  V/cm 

min 

V  =  5  V 

t  »  5.71x10  ^  s  drift  limited 

gate 

t'  ■  5.52x10  ®  s  diffusion  limited 
gate 

For  narrow  gate,  such  as  L=1  pm,  drift  velocity  saturation  sets  in  so 

3 

that  the  gate  transfer  time  decreases  with  L  rather  than  L  .  Thus, 

Equation  (7.9)  must  be  used  to  estimate  the  transfer  time.  The  results 

are  given  below. 

E  .  =  25  KV/cm 

min 

V  =  5  V 

t  ■  2.5x10  11  s  drift  limited 

gate 

t'  =  5.5x10  ^  s  diffusion  limited 
gate 

In  both  of  these  case,  the  E  .  under  the  gate  is  sufficiently  high  that 
the  charge  transfer  time  under  the  gate  is  limited  by  drift  rather  than 
by  diffusion. 

The  total  charge  transfer  time  is  the  sum  of  the  gate  and  gap 
transfer  times.  However,  for  overlapping  gates,  the  gap  transfer  time 
is  quite  small  compared  with  the  gate  transfer  time  and  can  be  neglected 
as  a  first  estimate. 

The  gate  transfer  time  given  by  Equation  (7.8)  shows  that  the 
ultimate  limit  is  affected  by  the  maximum  tolerable  electric  field  in 

gg  •  The  value  we  used  for  the  numerical  example,  lO"1  V/cm, 


the  silicon,  E 
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is  a  rather  high  field  which  will  give  i  electron-hole  pair  production 

rate  of  about  100  era  ^  per  electron  or  an  increase  of  electron  density 

of  0.01  in  a  gate  of  10  14  cm  length.  This  1%  increase  will  become 

significant  after  passing  through  100  charge  transfer  gates.  In 

addition,  the  noise  associated  with  the  random  pair  production  will 

also  limited  dynamic  range.  Thus,  a  smaller  value  of  E  :>r  applied 

3S 

voltage  between  adjacent  gates  may  have  to  be  used  which  would  reduce  the 
frequency  response  or  increase  the  transfer  time  of  the  device. 

Reliable  interband  impact  ionization  rate  data  at  low  electric 
field  are  not  available.  Some  measurements  have  been  obtained  at  low 
fields  using  junction  gate  field  effect  transistors  [48],  Charge 
coupled  devices  offer  alternative  structures  to  get  the  low  field  data. 

A  carousel  or  racetrack  three-phase  structure  has  been  designed  and 
fabricated  with  two-micron  gap  and  relatively  large  gate  length. 

The  top  view  of  a  device  is  shown  in  Figure  7.1.  The  circulating 
charge  packets  will  generate  electron-hole  pairs.  The  generated  holes 
will  flow  into  the  p-tvpe  silicon  substrate  and  produce  a  substrate 
current  from  which  the  hole  ionization  rate  can  be  obtained.  Since 
there  is  essentially  no  substrate  current  without  electron-hole  pair 
production,  the  substrate  current  measurement  can  detect  very  low 
hole  currents  at  low  impact  ionization  rates  and  low  electric  field. 

In  addition,  the  CCD  structure  gives  the  impact  ionization  rate  at  the 
Si-SiO,,  interface  where  additional  surface  scattering  mechanisms  may 
influence  the  electric  field  dependence  of  the  ionization  rate.  This 


project  has  just  begun  and  will  continue  in  a  new  contract. 
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